Impurity scattering in unconventional density waves 



Balazs Dora 

Department of Physics, Technical University of Budapest, H-1521 Budapest, Hungary 

Attila Virosztek 

- - - Department of Physics, Technical University of Budapest, H-1521 Budapest, Hungary and 

g : Research InstMute for SoUd State Phys.cs and Opt.cs, P. O.Bo. 49, H-1525 Budapest, Hungary 

o 

(N 



o 



o 



I 



X 



Kazumi Maki 
Max Planck Institute for the Physics of Complex Systems, 
Nothnitzer Str. 38, D-01187, Dresden, Germany 
(Dated: March 11, 2009) 



We have investigated the effect of nonmagnetic impurities on the quasi-one-dimensional uncon- 
ventional density wave (UDW) ground state. The thermodynamics were found to be close to those 
(N ■ of a d-wave superconductor in the Born limit. Four different optical conductivity curves were found 

depending on the direction of the applied electric field and on the wavevector dependence of the 
gap. 
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^ ; I. INTRODUCTION 



I , Recently a number of papers has been published investigating the different properties of unconventional density 



waves (UDW) under various conditions. The common feature of these systems is the zero average of the gap on the 
Fermi surface, resulting in the absence of any periodic modulation of the charge or spin density. Clearly this property 



^ makes UDW a very likely candidate for those systems in which clear thermodynamic signals of a phase transition are 
I— ' , detected without any obvious order parametero. From this the notion " hidden-order" follows naturally. 

Unconventional density wave formation is possible in a large variety of systems. In the quasi-one-dimensional case, 
^ ' which is the natural occurence of densitv^wavesH, we have investigated the basic propedjies of unconventional spin a: 
. charge density waves (USDW, UCDW)B and the related threshold electric field withcl and without magnetic ficl 
■ UCDW turned out to be relevant in the explanation of response of low temperature phase of quasi-one-dimensional 
, a-(BEDT,-TTF)2KHg(SCN)4 salts. In two dimensional systems, the different unconventional phases were elaborated 
' by OzakiQ. Among them, the d-density wave scenario which is a special case of UDW |(prbital antiferromagnet), was 
] proposed recently to describe the famous pseudogap phase of high Tc superconductoraJ. Since the original proposal, 
' several works have been published in which the properties of d-density waves were studied with the aim of testing 
the validity of the model by comparing to experimental data (see Ref. | rSfid the references therein). Also the ground 
state of certain heavy fermion materials were suspected to be USDWl1E£| which would simply explain the unsolved 
problem of micromagnetism. In the presence of magnetic field, the orbital antiferromagneillil and the spin nematic 
statellj were discussed as well in two dimensions. In three dimensional systems, the pseudogap phase of the transition 
'"^ metal oxides have atixacted significant attention and the staggered flux state was mentioned in the context of the 
^ , possible explanationgl3. p 

O ■ In this paper we extend our earlier analysiaJ on pure unconventional density waves to the presence of nonmagnetic 
impurities. Impurities are treated in the Born scattering limit since it works very well for conventional DW. Since 
the Fermi surface of quasi-one-dimensional systems mainly consists of two separate sheets, two different scattering 
processes should be taken into account: forward and backward scattering during which an electron remains on the 
^ ^ , same or moves to the other Fermi sheet, respectively. The thermodynamics are found to be similar to those of a d-wave 
' superconductor in the Born limit. Among the transport properties the quasiparticle part of the optical conductivity 
is evaluated. In the chain direction the phason couples strongly to the electromagnetic field, giving rise to massive 
collective modes in this direction. On the other hand, for electric fields applied perpendicular to the conducting chain, 
the conductivity shows only Fermi liquid renormalization, and our description is valid under these circumstances. 
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II. FORMALISM 



To start with, we consider the Hamiltonian of interacting electrons: 

^ = H ^(^)"^.<T°'k,^ + ^ l^(k,k',q)a+_^q^^ak,^aj,_q^^,ak',a', (1) 

^''^ k,k',q 
a, a' 

where aj^^ and ak,(j are, respectively, the creation and annihilation operators of an electron of momentum k and spin 
a. V is the volume of the sample. Our system is based on an orthogonal lattice, with lattice constants a, 6, c toward 
direction x, y, z. The system is anisotropic, the quasi-one-dimensional direction is the x axis. The kinetic-energy 
spectrum of the Hamiltonian is: 

^(k) = —2ta cos(fcj,a) — 2th cos{kyb) — 2tc cos{kzc) — fi. (2) 

In the second term of Eq. (|^) we consider the interaction between on site and nearest neighbor electrons as in Ref.i. By 
moving from Bloch space to Wannier space, the Wannier function is well localized, leading to a significant dependence 
of the interaction matrix element on the incoming electron momenta k and k'. Its antisymmetrized (therefore spin 
dependent) versionEJ is given by 

N ~ . ^ 

— F(k, k', q, cr, cr') = {U + 2_^{2Vi cosqiSi + 2 J, cos(fc, - + q,)5i + 

i 

-h2i?e(Fie^('=''+'=')*0 + 2Re{a{e'''^^^ + e'<^^ + e^C'^-J')'^' + e*(fe.+9.)'5.)))) + 

+5a.a' '^iVi - Ji){cosqi5i - cos{ki - kl + qi)Si), (3) 

■i 

where i — x,y, z and Si = a, 6, c, the different matrix elements involve the on site (U), nearest neighbour direct (V^), 
exchange (Ji), pairJnopping (Fi) and bond-charge (Ci) terms. This interaction is able to support a variety of low 
temperaturft_pliasea3, but we are only interested in unconventional DW (whose gap depends on the perpendicular 
momentum)B'li3. The latter can be either UCDW or USDW depending on the strength of the exdiapge and pair- 
hopping integrals. The single-particle electron thermal Green's function using Nambu's notation istjO 

G^ik,iu;n) = - dr(^,^',(k,r)^+(k,0))ffe'""^ (4) 
Jo 

where the Green's function is chosen to be diagonal in spin indices and the momentum space is divided into k and 
k — Q spaces (left- and right-going electrons) by introducing the spinors: 

ujn is the Matsubara frequency, Q = {2kp,Tr/b,7T/c) is the best nesting vector. The inverse of the above Green's 
function is obtained as 

G^^(k,iw„) = iujri - ^(k)p3 - Acr(k)pi, (6) 
where pi (i = 1,2,3) are the Pauli matrices acting on momentum space, Ao-(k) satisfies the self-consistent equation: 

W = 7 E nk',k,Q,a,a')(at,,,ak'+Q,.). (7) 

k',cr' 

In order to describe USDW, we assume A as an odd function of the spin (A^- = — A_(j). Assuming to be an even 
function of the spin, we would have UCDW. From now on, we will drop the spin indices since they are irrelevant for 
most of our discussion and most of our results applies to both unconventional charge and spin density waves. The 
spin indices will be reinserted wherever necessary. With this, the gap equation reads as 
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where E(k) = •^/^(k)^ + A(k)p, A(k) = Acr(k) and the kernel of the integral equation is diagonal on the basis of 
the leading harmonics asL 

-P(k 1) Pq Pi P2 -P3 P4 

— = ^ + ^ cos{kyb) cos{lyb) + sin(fcy6) sm{lyb) + 003(^^0) cos(Z^c) + sin(/c^c) sin(;zc). (9) 

The Pi coefficients are linear combinations of the interaction matrix elements. As a consequence of the general form 
of the kernel, the gap will be of the form 

A(l) = Ao + Ai cos{lyb) + A2 sm{lyb) + A3 cosH^c) + A4 sin(Z^c). (10) 

From now on we assume that only one kind of gap among the five possible candidates, whose transition temperature 
is the highest, opens and persist all the way down to zero temperature. For example we find that USDW is stable wi 
respect to UCDW if JyTFy > 0, where the upper (lower) sign refers to a ky dependent gap function of cosine (sine)l 
The thermodynamic and transport properties of such a system has been worked out in Rcf. ^. In the foUowings 
we shall discuss the effect of impurities on UDW and determine the behaviour of the basic physical quantities. The 
interaction of the electrons with nonmagnetic impurities is described by the Hamiltonian: 

^1 = -^ E e-^^^^^+{^ + ^)U{Il,)^a{\^), (11) 



V, 



_ / UiO) C/(Q)e-Q^. \ 



(12) 



Rj is the position of the j-th impurity atom. The explicit wavevector dependence of the matrix elementsEH13 is 
neglected since no important changes are expected from it. The usual method of treating the impurities is to average 
over their position in real space, and step into the wavevector space afterwardsllZlo. Instead, we follow a rather 
unorthodox way: working in the Fourier space and averaging when needed. It is clear from the exponential prefactor 
in Hi, that only diagrams containing impurity scattering with momentum conservation at each impurity atom have 
finite expectation value after averaging over the position of the impurities, and translational invariance is regained. 
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FIG. 1: Self energy corrections due to impurity scattering. The solid line denotes the electron while the dashed line is for the 
electron-impurity interaction. Dashed lines coming from the same cross represent successive scattering of the electron on the 
same impurity. 

This method can be extended into any order of impurity scattering as we will d ernonst rate it in the foUowings. As to 
the diagrams, we will take into account only non-crossing, ladder type diagramslljllllc3. Recently the applicability of 
this approximation in 2 dimensions has been questioned apd a new method has been invented in order to consider all 
types of diagrams (i.e. those with crossing impurity lincs)Eil. As a result, a novel type of ni/|2cj|(ln^ I'^/^l + ('''/2)^) 
additional density of states was found, which could not be obtained by the non-crossing apprcj ^dm ation, Ui is the 
impurity concentration. But for 1 and 3 dimensional systems the usual technique looks sufficientEHlcJ. To start with, 
we will evaluate the self energy corrections caused by Eq. (|l^) at every order. This can be visualized in Fig. |l] and 
is given by 

ER(k,io.„) = SR(ic^„) = (u{-R) + C/(R) J 0G(p, io.„);7(R) + 

+C/(R) J 0G(p, iiv^WiR) J ^G(p', zc^„)C/(R) + . . = n,f/(R) + [/(R) J 0G(p, ic^„)ER(zL^„), (13) 

where the self energy correction turns out to be momentum independent and the R index in SR(iw„) means the 
position of an impurity over which the average will be taken in the foUowings. Eq. ( p^ can be solved easily, and the 
result is: 

ye ,f_Ui-g_U2e-''^^ + f\ ni 

2.R^*a;„j l^C/^e'QR + j Ui - g / ([/i - g)2 - |/|2 _ jc/^p _ ([/Je-Q^ + Us/e^Qi^) ' ^ ^ 
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where Ui = U{0)/{U{idf - \U{q)\^) and U2 = U{q)/{U{Qf - \U{CtT) and 



'^'^.G(p,zc.„)= f 4 { ). (15) 



Expanding ER(ia;„) in powers of the exponential terms in the denominator of Eq. (^^; the space average can be 
performed and the self energy matrix is obtained as 

i](^.„) = fg!H SSr-n (16) 



T.'i{iu!n) Si(zu;„) 



and its matrix elements are given by 



ijiiiujn) — ni , , 17 

E2(zc.„) = ^ = ( U{Of - |C/(Q)|2 + (^8) 

v/i?2-4|;7(Q)/|2 V^2_4|[/(Q)J-|2y 

S3(«w„)/ = S2(iu;„)/, (19) 

where D = 1 — 2(7J7(0) + (5^ — (0)^ — |C/(Q)p). This result is valid only for a certain range of parameters due 

to the expansion. On the other hand, one can deduce an expression involving the different matrix elements of the self 
energy where the average can be performed rigorously: 

Ei(l - 2C/(0).g + {g^ + \m{U{Qf - |C/(Q)n) - + S3/ + n,)(C/(0) - g{U{Qf - UiQf) = 0, (20) 

and this equation is satisfied with the previously obtained Si, S2 and S3 even outside of the validity range of the 
expansion. Of course this cannot be regarded as a proof but we can trust in the usefulness of this calculation outside 
the validity range. Moreover in a normal metal this result gives back the known resultEfl. These formulas apply also 
to superconductor|S_mith.pmiaor change (C/(R) = U{0)pi), and the self energies in the Born and unitary limit are 
obtained correctlyE£lE5BEat3. Wcrireat our UDW system in the Born scattering limit since conventional DWs are 
commonly investigated in this limited. The interaction gives rise to the self energy, which is in the Born-approximation 
(considering only the lowest order terms): 



S(k, toj,,) = y E ^ E U{Il)G{k - q, tu;n)U{R), (21) 

q R 

where the summation is the only remaining operation from averaging over the impurity atoms. From this, one obtains 
for a DW 

Gi^,,^^^)^-'.^^^l±i^}Pl±h^, (22) 

where both the frequency and the gap are renormalized in the conventional case: 

ri + Fa 

w„ = w„ (23) 

A = A„ + ^ , ^" . (24) 
■ + A2 



Fi — 7rni|C/(0)p.g(0) is the forward scattering, F2 = 7rni|[/(Q)pg(0) is the backward scattering parameter. Ui is-jthe 
impurity concentration, 17(0) is the density of states per spin in the metallic state. As in other similar problemst3, it 
is convenient to introduce the quantity u„ — a)„/A„, which relates to physical quantities: 

u;„ = Au„ I 1 - a ] ) , (25) 



5 



r = Fi + -ij-, a = r/A is the pair-breaking parameter. As opposed to this, in unconventional DW self energy 
corrections from impurities does not renormalize the gap, only the Matsubara frequency: 

Ti + r2 a)„ A \ 

= W„, A 




IT ^Col + A2 \^Col + ^^) 
A„(k) = A(k) = Asin(6fcj^) or Acos(&fcy). (26) 

This is written in a more useful dimcnsionless form: 

T ^, / 1 \ \ 

(27) 

where T = (Fi + F2)/2, a = F/A, m„ = uin/^ and Fi and F2 are the same quantities as in a conventional DW, K{z) 
is the complete elliptic integral of the first kind. Here the combination of the scattering rates is different from the 
conventional DW's case due to the lack of renormalization of the order parameter. We choose the Born scattering 
limit because this limit works very well for conventional DW. We believe that by neglecting the explicit wavevector 
dependence of the impurity matrix elements, we made a useful approximation as far as the character of the physics is 
concerned and we are able to capture the characteristic changes caused by inap.u£j|tira_ However in order to describe 
very fine, characteristic phenomena to DW such as the threshold electric fieldc§Ejl2j'Lil, we cannot use simple s-wave 
scatterers as it is shown in Refs. ^,|[ 



III. THERMODYNAMICS OF IMPURE UDW 

Since the thermodynamic properties of a pure UDW are identical to those of a d-wave superconductoitJty and 
the impurity effects on a conventional DW are similar to those in s-wave superconductors, we expect very similar 
behaviours to those in a d-wave superconductor treated in the Born limit. However, the main difference is that we 
distinguish two different scattering processes (forward and backward scattering) while in the superconducting world 
there is only one. Consequently the different combinations of the F's are far from being trivial. The gap equation is 
obtained as 

' - ? (to) ^ - " (™) A) • 

where E{z) is the complete elliptic integral of the second kind. The change in the transition temperature is given by 
the Abrikosov-Gor'kov formula: 

-ln(^]^^(l+p)-^j(l], (29) 



TcoJ V2 J \2, 

where Tc and Tc^ are the transition temperature of the impure and clean system, respectively, p = T/2t:Tc, ^{z) is the 
digamma function. Note that this formula is also valid for any kind o£-unconventional superconductor in the presence 
of impurities considered either in Born or in resonant scattering limitS. The critical impurity scattering rate is given 
by 

The gap maximum is the same as the one of a d-wave SC in the Born limitS: 

1"a^ = ^X r{K-E)(E^K-^)dx + 2{sin^r^rs\.-^), (31) 
A(0,F) TT-^AJco \ l + x-'J £--"• 



smy 



where (...) means ^ f^^ dy . . . , the argument of K and E reads as . Cq is the value of u„ at zero frequency 



2 



C2 = -aK , , (32) 
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vanishing as the impurity scattering parameter disappears hke Co = 4exp(— 7r/2a), while for large a: Cq = a as in 
Ref. pq. Close to T^, A vanishes in a square-root manner as does usually in mean field treatments: 



= m-T^r^TTT^^^^A^ fl - - ) • (33) 



2 P 3 

Close to absolute zero, the following formula is obtained: 

A(r) = A(0)-y^ f--l) (1-3 (— / {K-E)x 



Co 

2 X , / ^2 XXX -1 



^ {E-KYf^)d- + Co^^ + C^o[E^K^^]]] T\ (34) 



where the decrease of the pure case turned into a faster one. 
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FIG. 2: A(0,r)/Aoo (dashed line), Tc/Tcq (solid line) and iV(0,r)/7Vo (dashed-dotted line) are shown as a function of V/Vc 
for an unconventional density wave. 



From this, one can assume that the effect of the impurity scattering in the limit of low temperatures is to reduce 
the power-law exponent by one. As a-^esult we expect the exponent of temperature to be the same as those in a 
conventional DW in the gapless regionll3'Ej. The analogy looks obvious since in neither of these two systems there is 
a lower bound of the excitation energy. The correspondence works only at low temperatures T <C Tc when the only 
energy scale is the temperature. Now we derive expressions for the grand canonical potential and for the specific iieat. 
In doing this, we use the well-known relation involving an integral over the coupling constant of the interactiouEj: 

/■I 

f]-0o = —{XH,„t), (35) 

where Hint is the interaction causing the phase transition. This formula gives us the thermodynamic potential 
difference between the normal and the DW phase. Since we work on a grand canonical ensemble, the appropriate 
thermodynamic potential at T = is obtained as: 

m = -NpiO) - ^A^CoV^+ y + ^ - ^ - ^) - ^TT^) ■ (36) 

At small r, the leading correction is the last integral, enhancing the potential as in the normal SDW case. The low 
temperature specific heat reads as 

cm = ^5(0)^T. (37) 
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This expression also reaches the normal state value with increasing F. The specific heat jump is: 



2 P 3 

In Fig. ^ we show A(0, F) and Tc as a function of the scattering rate. 



(38) 



IV. DENSITY OF STATES IN UDW 




FIG. 3: Density of states plotted as a function of the reduced energy for different scattering amplitudes: a = 0, 0.01, 0.05, 0.1, 
0.5 and 1 with peakposition at u = A from top to bottom. 



By use of the Green's function, the density of states per spin is given by: 

= -;-^5]lmTr(G«(k,c.)) =g(0)ilm(u), 

k 

where u = iun{iujn = iu + iS). After some algebra, the low energy behaviour reads as: 



iV(L^) = 5(0)^ ( 1 



8E^ 



f<2 



f 



(39) 



(40) 



The residual density of states (i.e. the DOS at the Fermi energy) is finite at any finite F, disappears exponentially as 
F goes to zero, but takes the normal state value as F approaches to infinity. Since A^(0) is almost zero for F < O.SFf,, 
we do not expect relevant changes in the static quantities (such as the specific heat, the spin susceptibility at T ^ 0) 
at low impurity concentrations. The notion "gapless" makes no sense in this case since even in pure UDW the gap 
vanishes at the Fermi energy leading to the possibility of arbitrary small energy excitations. At the value of the order 
parameter, the divergent peak of the pure system is broadened arid.iV(a;) is always finite as a result of the impurities, 
shown in Fig. ^. Compared to the DOS of the conventional DWEj, the states below the gap maximum are filled in, 
and the peak at A is disappears more rapidly as a increases than in the case of momentum independent gap. At high 
energies it reaches the normal state value as: 



N{lo) = .9(0) 1 + 



A2 a;2 



4 (w2 + F2)2 

In Fig. ^, we show the F dependence of the residual density of states. 



(41) 
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V. DENSITY CORRELATOR 

We turn our attention to the behaviour of the static, long wavelength density correlation functionll using the the 
thermal Green's function: 

^o(^) = -^ E Tr(G(p,k,zu;„)G(k,p,ic^„)) (42) 

p,k,CT,n 

where the overline means averaging over the position of the impurity atoms. This requires calculating the averaged 
Green's function and the vertex corrections, since the average of the product of two Green's functions is not equal to 
the product of the averaged Green's functions. In the following we focus on the vertex corrections, A(p,za;„). With 
this, our equation becomes simpler: 

^o(r) = E Tr(G(p,za;„)A(p,za;„)G(p,z^„)). (43) 

p,cr,n 

In the standard ladder approximation the vertex corrections are determined by the integral equation: 

A(p, za;„) = 1 + y E ^ E WG(q, *'^n)A(q, zc^„)G(q, iu;„)U{R), (44) 

q R 

which is shown in diagrammatic language in Fig. ^. Assuming A(p, zcij„) = A{iujn), and making the following ansatz: 




FIG. 4; The vertex correction in the Born limit is shown. The dot is the vertex function, the filled triangle represents the 
vertex correction due to impurity scattering. 



the vertex corrections can be obtained: 



Aa = 



Ai(iw„) A2{iLu„ 
A2(iw„) Ai(iw„) / ' 



2 K - E 
1 a- 



Substituting this to Eq. (p3|), the susceptibility reads as: 

/ 

xo(r) = 2<?(o) 



K-E 



^ 2 ^ ^ul + \ 



\ 



A/3. 



v 



1 - a- 



K- E 



At zero temperature it equals to the total density of states at the Fermi surface 



Xo(0) = 2g(0)- 



F ■ 



In the low temperature limit we obtain: 



T 



Xo(T) = 2,(0)^ A(T) + ^(f-l + ^^^^ 



(45) 

(46) 
(47) 



(48) 



(49) 



(50) 



9 




Close to Tc, a similar expression to the normal SDW describes the susceptibility: 

\ 2 P 3 

We refrain from the evaluation of the Q-th Fourier component of the density correlation function because in UDW 
this is not the quantity which signals the phase transition, RPA corrections, will not lead to divergence, since the 
dominant unconventional channel does not couple to charge or spin densityrl. In conventional CDW or SDW, the 
Q-th Fourier component of the charge density or the spin density turned out to be the order parameter of the phase 
transition, respectively. As opposed to this, in the unconventional scenario, the following phases and related order 
parameters are found: 



phase 


gap 


order parameter: the Q-th 
Fourier component of the 


UCDW 


A cos{bky) 


electric current density 


UCDW 


A 8m{bky) 


kinetic energy density 


USDW 


A cos{bky) 


spin current density 


USDW 


A sin(fefcj,) 


spin kinetic energy density 



These phases are akeady known as orbital antifcrromagnetO, bond-order wavcQ, spin nematic stateEHl and axial 
spin bond-order waveO, respectively in the context of the two. dimensional Hubbard model. Generally these order 
parameters can be called as the effective charge or spin densityEl. The autocorrelation function of the above quantities 
will be divergent at Tc in the corresponding phase, because these are the relevant quantities from the phase transition's 
point. 



VI. OPTICAL CONDUCTIVITY 



The optical conductivity contains relevant informations about the possible excitation of a system. Since in real 
materials impurities are always present, the evaluation of the optical conductivity in impure systems is of prime 
importance. As it is known, the electriGal conductivity of a conventional DW is divided into a pair-breaking (interband) 
and a normal (intraband) contributionE^. Hence a Lorentzian like normal contribution appears at all the frequencies, 
while the pair-breaking term is zero as long as a; < 2A. This separation can be done in the unconventional case, 
although here both processes contribute to all frequencies due to the finite density of states at the Fermi energy. 
Introducing two notations: 

/■°° / Blx + uj) 0x\ 

In{uj)^J ftanh--— ^ -ta.nh-^ \Re{F{u{uj + x),u{x)) - F{u{lj + x),u{x)))dx, (52) 

Ipb{Lu)^[ tanh ^'^^ ~ Re{F{u{oj - x),-^il{x)) - F{u{uj - x), -u{x)))dx, (53) 
Jo 2 

the conductivity is given by: 

Reaaa = -e 9{0)v^- ^ , (54) 

Att lu 

where Vx = vp, Vy = ^/2btb and Vz = V^ctc. The different F{u, u') hmctions and the dc conductivities are commented 
in the foUowings: 

i. A(k) — Acos(fcy6), a~y: 



F{u,u') = 
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This is the simplest case, the vertex corrections vanish due to the mismatch of wavevector dependence of the velocity 
and the gap. As the scattering strength enhances, it becomes the dominant energy scale and the curves take more 
and more the form of a Lorentzian as it can readily be checked in Fig. |^. The DC conductivity is calculated at T = 0: 

oty- - e^mvl^ [e^^I - ^) . (56) 




FIG. 5: Real part of the electric conductivity in the y direction for A(k) = Acos(6fcy) is plotted as a function of the reduced 
energy for different scattering amplitudes: q = (dotted line), 0.1 (solid line), 0.5 (dashed line) and 1 (dashed-dotted line), 
Fi = 2F2. 



ii. ^(k) = Asin(fcj,6), a — y: 

1 



F{u,u') 



+ 



ri 1 



, 2 u' 



K [ -uu' + 3 + y 

,2 



+ 



-.K' - 



u 



-.K 



Att (u + u') 



A2 



Ti 1 / u 



(57) 



-.K 



the third row of the equation comes from the vertex corrections. As F increases, the peak at 2A is broadened and 
moves closer to zero frequency. The DC conductivity is obtained at T = 0: 



yy 



Cl{K-E) 



^^:^JCl + l + Tl{K-Ey 



(58) 



where the second term in the denominator is clearly the effect of the vertex corrections. The conductivity is shown 
in Fig. |. 



iii. A(k) = Asin(fcj^6) or Acos(fcj,6), a = z: 

1 



F{u,u') 



2{u'' - w2) 



2^/l -u^E - 2\J\ - u'^E' + K 



.u'(u~u') ^^u{u~u') 
' ) ' + K \ ' - 



(59) 
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a;/Ao() 

FIG. 6: Real part of the electric conductivity in the y direction for A(k) = Asin(6fcj,) is plotted as a function of the reduced 
energy for different scattering amplitudes: a = (dotted line), 0.1 (solid line), 0.5 (dashed line) and 1 (dashed-dotted line), 
Fi = 2r2. 



the vertex corrections vanish because the velocity depends on different perpendicular wavevector component (fc^) than 
the gap (ky). As F increases, the peak at 2 A is broadened and moves closer to zero frequency. The DC conductivity 
is obtained at T = as 



dc 



2e'g{0)vl 



E 



(60) 



The optical conductivity is usually the same in the x and z direction apart from constant factors, since the velocity 
in these directions does not interfere with the gap. But in the presence of impurities this general relation does not 
hold any more due to the presence of different vertex corrections. A very similar breakdown of equality is found in 
the relation between the static spin susceptibilitv-and the condensate density (ps = 1 — Xo/Xn)i which are not related 
to each other if impurity scattering is considereclljcj. The conductivity is shown in Fig. 



CM 

b 




j/Ao 



FIG. 7: Real part of the electric conductivity in the z direction is plotted as a function of the reduced energy for different 
scattering amplitudes: a = (dotted line), 0.1 (solid line), 0.5 (dashed line) and 1 (dashed-dotted line), Fi = 2F2. 



For the sake of completeness we present the result for the quasiparticle part of the conductivity in the chain direction 
keeping in mind that collective modes also appear in this direction. 
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IV. 



A(k) = Asm{kyb) or Acos{Kyb), 



F{u,u') 



ttA 



2(ri -Ta) 



Att 



-if 



u{u' — u) ,u'{u' — u) 



^^^^^ 



K' 



2E'V1 



(61) 



This formula gives the quasiparticle part of the optical conductivity in the chain direction, although collective modes 
also show up here significantly modifying the conductivity. The consideration of impurity scattering anjd-|eallective 
modes (even in the simplest random-phase approximation) together is a very difhcult task to deal withEIrEa and is 
beyond the scope of the present investigation. The DC conductivity is obtained at T = 0: 



dc 



E 



ATrycf+l- {Ti-T2)E 



(62) 



The conductivity seems to transfer more and more spectral weight to the zero frequency peak with growing impurity 
scattering rate, transforming the curve into a Lorentzian like one (Fig. ^). 



CM 



b 

•A 




tj/Aoo 

FIG. 8: Real part of the electric conductivity in the chain direction is plotted as a function of the reduced energy for different 
scattering amplitudes: a = (dotted line), 0.1 (solid line), 0.5 (dashed line) and 1 (dashed-dotted line), Fi = 2F2. 

The dc conductivities are shown in Fig. ^ at T = as a function of the impurity scattering parameter. In the 
perpendicular direction, the dc conductivities take the same value at the critical scattering parameter, while the dc 
conductivity in the chain direction is exactly 3/2 times larger a s fol lows from Eq. ( |63|) and ( p^ ) in the uj = Q limit, if 
Fi = 2F2. In spite of the similar thermodynamics of d-wave SCat3 and UDW, the transport properties of these two 
systems are completely different due to the distinct coherence factors coming from the different condensates. In a SC, 
there is always a Dirac delta peak at zero frequency, and the stronger the impurity scattering, the larger the spectral 
weight of this peak transferred to the finite frequency part of the conductivity. In UDW, the Dirac delta contribution 
disappears as soon as any finite impurity concentration is present, and the areas under the different curves are equal, 
but their form approaches those in the normal metal as F enhances. 

The normal state electric conductivities are given by the usual Lorentzians: 



Recr;:,^(c^) = e^5(0)2w| 



2F2 



Recryy,^z(w) = e^g{Q)2vl 



W2 + (2F2)2^ 

2 ri + F2 



^c^2 + (ri + r2) 



(63) 
(64) 



In the chain direction only backscattering can cause current damping as it is known from transport theory, which is 
manifested in the absence of the forward scattering parameter in Yieaxx{^)- 
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FIG. 9: The dc conductivity plotted at T = as a function of the reduced scattering rate for Fi = 2r2 = 4r/3 for a cosinusoidal 
(sinusoidal) gap in the y direction: solid (dashed line), in the z direction: dashed-dotted line and in the x direction: dotted 
line. 

VII. CONCLUSION 

We have studied the effect of nonmagnetic impurities in unconventional density waves. In this respect there is no 
difference between USDW and UCDW due to the spin independence of the interaction with impurities. In s-wave 
superconductors nonmagnetic impurities have no influence on the thermodynamics of the system, while impure d-wave 
superconductors suffer important changes. This is known as Anderson's theorem, but equivalent conclusion has been 
reached independently by Abrikosov and Gor'kov. It says that if a static perturbation does not break the time-reversal 
symmetry and docs not c;ause a long-range spatial variation of the order parameter, the thermodynamic properties 
of the superconductor remain unchanged in the presence of perturbation. As opposed to this, any kind of DW is 
destroyed in the presence of impurities, although the identity of the thermodynamics of s-wave superconductor to 
conventional DW and d-wave superconductor to unconventional DW is well established without impurities. Impurities 
have a pair breaking effect on the condensate, resulting in a universal formula between the transition temperature and 
the scattering parameter, named after Abrikosov and Gor'kov. It seems to be valid for superconductors with all kinds 
of symmetries and now for density waves as well, independent of whether the Born or the resonant scattering limit 
is taken. Since conventional DW were studied in the Born limit, we found sufficient to use the same approximation 
for the unconventional scenario. We have examined the system with the standard non-crossing approximation, and 
calculated the self-energy corrections for infinite order in the scattering potential, but only the lowest non-trivial 
correction was retained for the Born limit. The thermodynamics of UDW were found to be very similar to the one in 
d-wave superconductors with nonmagnetic impurities, but the existence of two different types of scattering processes 
(forward and backward) was called for in the microscopic theory. In unconventional DW, at any finite scattering 
strength the valley of the density of states at the Fermi energy is filled in, leading to normal electron like behaviours 
very close to absolute zero, but the reduced density of states compared to the normal state bears the effect of the 
condensate. The order parameter does not get renormalized due to impurities because we assumed s-wave scattering 
for simplicity. The specific heat increases linearly with temperature due to the finite density of states at the Fermi 
energy. Interestingly, impure UDW was found to be very similar to the gapless region of conventional DW very close 
to the critical scattering rate as long as the temperature exponents are concerned close to absolute zero because of 
the absence of any finite lower barrier of the excitation energy. 

But at the transport properties all the similarities ended. The optical conductivity in the chain direction is domi- 
nated by the phason contribution, and incorporating the effect of impurities in the theory is beyond the scope of this 
study. Instead we concentrated on the perpendicular direction. In the optical conductivity, self energy and vertex 
corrections were taken into account in the ladder type non-crossing approximation. Depending on the symmetry of 
the order parameter and the chosen direction, four qualitatively different curves are deduced, although a^x is certainly 
dressed by collective modes due to coupling to the phason propagator. In the perpendicular directions, the possibility 
of low frequency excitations rapidly increases, transferring increasing amount of spectral weight to a; = 0. The dc 
conductivities at T = sharply differ from each other hence they can help to provide one with decisive conclusion 
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when comparing these results to experimental data. 
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